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The Collins fragmentation function: a simple model calculation 
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The Collins function belongs to the class of the so-called time-reversal odd fragmentation functions. 
Being chiral-odd as well, it can serve as an important tool to observe the nucleon's transversity 
distribution in semi-inclusive DIS. Due to the possible presence of final state interactions, this 
function can be non-zero, though this has never been demonstrated in an explicit model calculation. 
We use a simple pseudoscalar coupling between pions and quarks to model the fragmentation process 
and we show that the inclusion of one-loop corrections generates a non-vanishing Collins function, 
therefore giving support to its existence from the theoretical point of view. 
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Our understanding of hadronic physics depends 
strongly on what we know about the parton distribu- 
tion and fragmentation functions, which are universal, 
process-independent objects. While in the past several 
experiments provided us with considerable information 
on the parton distributions, our knowledge of the frag- 
mentation functions is still rather limited. 

A lot of attention has been devoted to the class of the 
so-called time-reversal odd (T-odd) fragmentation func- 
tions, among which the Collins function [Q,^ is the 
most prominent example. In contrast to a naive expecta- 
tion, these functions can be non-vanishing because time- 
reversal invariance does not impose any constraint [^]||] 
(in the literature they are sometimes referred to as naive 
T-odd functions Q). This is a consequence of possible 
final state interactions in the fragmentation process, giv- 
ing rise to a non-trivial phase (imaginary part) We 
remark that final state interactions should not be consid- 
ered exclusively as reinteractions of the outgoing hadron 
with the rest of the jet, as we will discuss in more detail 
later. 

The Collins function describes the fragmentation of a 
transversely polarized quark into an unpolarizcd hadron 
(e.g. a pion), i.e. the process q* — > hX. The introduction 
of this function requires taking into account the quark's 
transverse momentum. Although the Collins function 
in itself deserves attention towards understanding the 
physics governing the fragmentation process, it is partic- 
ularly relevant because in semi-inclusive DIS it can serve 
as the chiral-odd partner needed to access the transver- 
sity distribution of the nucleon, hi, which is chiral-odd 
as well. The transversity is a crucial property of the 
nucleon, carrying information on its spin structure com- 
plementary to what we can deduce from helicity distribu- 
tions. Because of its chiral-odd nature, hi is suppressed 
in inclusive DIS, and therefore it has remained essentially 
unexplored on the experimental side. In the case of one- 



particle inclusive DIS, where one detects only a single 
hadron in the final state, the Collins function provides 
the only possibility to probe hi at the leading twist level. 

Even if the Collins function can not be discarded in 
general on the basis of time-reversal invariance, it has 
been suspected to vanish since the effect arising from final 
state interactions may average out |^. Moreover, has 
been proven elusive to previous modeling attempts [Q. 
This suggested to turn the attention to the experimen- 
tally more challenging detection of two hadrons in a semi- 
inclusive measurement |6|j^, as an alternative method to 
measure the transversity. Until now, no ah initio calcu- 
lation has ever displayed the possibility of generating a 
non-zero Collins function in the framework of a simple, 
time-reversal invariant model. 

On the experimental side, the HERMES collaboration 
reported the first observation of a single-spin asymmetry 
in semi-inclusive pion production 0, which could be in- 
terpreted as arising from the contribution of the Collins 
function Analogously, large single-spin asymmetries 
in the process p^p — > ttX |^ could also be explained by 
means of the Collins function ||l^]. The situation is far 
from clear, but further investigations are among the pri- 
orities of the HERMES 0, COMPASS (l| and possibly 
eRHIC programs. 

In these circumstances, producing a non-zero Collins 
function in a simple yet consistent model is an interest- 
ing and relevant result, which we are going to present in 
this letter. To this end we describe the fragmentation 
process by a pseudoscalar coupling between quarks and 
pions in the one-loop approximation. A similar calcula- 
tion has been suggested earlier |l3[ , but never carried out 
explicitly. For what concerns possible phenomenological 
applications, we do not pretend our calculation to be re- 
alistic, but we think that it can shed light on the identity 
of T-odd fragmentation functions and can conclusively 
affirm that there are no theoretical reasons to believe the 
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Collins function to vanish. 

Considering the fragmentation process q*{k) — > 
t:{p)X, we define the Collins function, which depends on 
the longitudinal momentum fraction z of the pion and 
the transverse momentum fc^ of the quark, as ^ 
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dfc+ Tr[A(fc,p)icr'-75] 



FIG. 1. Lowest-order unitarity diagram describing the 
fragmentation of a quark into a pion. 
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dfc+ Tr[A(fc,p)7" 



with TO^ denoting the pion mass and e^? = e'^ + (we in our model at tree level becomes 
specify the plus and minus lightcone components of a 
generic 4- vector according to = (a"±a'^)/'\/2). The 
correlation function A{k,p) in Eq. (|l|), 
links, takes the form 
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omitting gauge 



Di{z,k'^) = - 



Z 16tT^ (fcj -(_ 77j2 



(5) 



(6) 



A{k,p)^Y. 

X 



|V(e)k,^)(7r,X|^(0)|0). 

(2) 



To describe the matrix elements in the correlation 
function, we use a pseudoscalar coupling between quarks 
and pious given by the interaction Lagrangian 



Ci{x) = ig q{x)-i^q{x)TT{x) , 



(3) 



which is in the spirit of the Manohar-Georgi model [|5| . 
This is clearly an oversimplified approach to the frag- 
mentation q* — > TiX but the model contains the essential 
elements required for our discussion. In particular, it 
is time-reversal invariant. One could also perform the 
calculation in a chirally invariant model by including a 
scalar a field as well as taking quark fiavors properly into 
account. In our view, this is certainly necessary to im- 
prove the result obtained here for a better description of 
the phenomenology, but the essential result of a non-zero 
is already evident without going to the complica- 
tions arising from the inclusion of the a particle and the 
discussion of various fiavors. 

Using the Lagrangian in Eq. (||), at tree level the frag- 
mentation of a quark is modeled through the process 
q* —f Trq. The corresponding correlation function can 
be represented by the unitarity diagram in Fig. |l] and 
reads explicitly 



^{o){k,p) = - 
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where m represents the mass of a constituent quark. Us- 
ing this correlation function, the unpolarized fragmenta- 
tion function 



*Note that this definition of slightly differs from the 
original one given by Collins B. 



In the case of TOtt = we recover the result already ob- 
tained by Collins Contrary to Di, the Collins func- 
tion at tree level is zero. This is not surprising because 
at tree level no final state interaction appears in the frag- 
mentation process, which is supposed to be the origin of 
the T-odd fragmentation functions like H^. The situa- 
tion changes when we proceed to one-loop corrections, as 
we explicitly show in the following. 

A consistent one-loop calculation of the fragmentation 
process requires the evaluation of all the diagrams shown 
in Fig. 2. In the actual calculation only the asymmetric 
diagrams (a-d) contribute. The contributions from the 
diagrams (e), (g) and (h) vanish individually when pro- 
jecting out the Collins function, while the contributions 
of diagram (f) and its hermitian conjugate cancel each 
other. 

It should be noted that not all the asymmetric 
(interference-type) diagrams contribute. For example, 
the diagram (f) and its hermitian conjugate represent 
the interference between two different tree level ampli- 
tudes describing the process q* tttt^. However, their 
contributions to the Collins function cancel each other 
since the involved interfering amplitudes are purely real. 
In calculating Hi , such a cancellation between contribu- 
tions coming from a diagram and its hermitian conjugate 
takes place whenever the involved amplitudes are purely 
real and is a model independent feature. 

In our computation, the relevant components to be 
included are the self-energy correction for diagrams (a) 
and (b), and the vertex correction for diagrams (c) and 
(d) . These corrections are sketched in Fig. || and can be 
analytically expressed as 
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[{l - kY - m^][l^ - ml] ■ 
The functions S(A;) and T(k,p) can be parametrized as 
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FIG. 3. One-loop self-energy and vertex corrections. 
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FIG. 2. One-loop corrections to the fragmentation of a 
quark into a pion. 



I](fc) = Af^Bm, 
T{k,p) ^C + D^ + Eft + F^^. 
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The real parts of the functions A, B, C etc. are UV- 
divergent and require in principle a proper renormaliza- 
tion. Though our model is renormalizable, we do not 
have to deal with this question at all, since only the 
imaginary parts of the loop diagrams will turn out to 
be important. 

The contributions to the correlation function generated 
by the diagrams (a) and (c) are given by: 
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The contributions from diagrams (b) and (d) follow from 
the hermiticity condition: A|J^^(fc,p) = 7°A|°j^(fc,p)7'^, 

A|^;(A:,p) = 7°Agt(A:,p)7°. 



Summing the contributions of the four diagrams and 
inserting the resulting correlation function in Eq. (|l|), we 
obtain the result 
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Thus the actual value of the Collins function in this 
model depends only on the imaginary parts of the co- 
efficients defined in Eqs. (|-0). The lack of an imag- 
inary component in these coefficients would inevitably 
result in a vanishing Collins function. We can compute 
the imaginary parts by applying the Cutkosky rule to 
the self-energy and vertex diagram of Fig. H. In this way, 
as mentioned before, we can avoid the issues related to 
renormalization, which affect only the real parts of the 
diagrams. Explicit calculation leads to 
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(14) 



Im(_D-|-£;-|-TOF) = --^TO,- 

Stt^ A(k2,to2,to^) 

X [h + (fc' - to2 - 2to2 )/2] , (15) 
where we have introduced the so-called Kallen function. 



X{k^,m'^, ml) 
the factors 



[fc2 — (to -I- TOtt) — (to, - 



dH5{P -ml)5{{k~lf -IT?) 



-^^A(fc2, to2, to2 ) e{k'^ - (to + m^)^) 
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and 
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m^) 5{{k - If - m^) 



{k — p — If — m? 



2^A(fc2, m? , TO,2 ) 
X e{k^ - (to + m^Y). 



A(/s2, to2, to2) 



k'^m? — {m? 



to2)2 
(17) 



These integrals are finite and vanish below the thresh- 
old of quark-pion production, where the self-energy and 
vertex diagrams do not possess any imaginary part. 
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FIG. 4. Interpretation of the one-loop corrections as final 
state interaction. 



Thus Eq. (|lj) in combination with Eqs. (|14|-|17D gives 
an expUcit non-vanishing result for the Collins function. 

We point out that Collins suggested the idea of dress- 
ing the quark propagator as a possible mechanism to 
produce a non-zero Here, we have not only sup- 

ported this conjecture by means of an explicit one-loop 
calculation, but we have also shown that the contribu- 
tions due to the self-energy and the vertex corrections do 
not cancel each other. 

To avoid possible confusions, we would like to make a 
few remarks on the issue of the final state interaction. 
One-loop corrections in our model can be viewed as con- 
taining a specific example of final state interaction, where 
the pion, after being emitted and before being detected, 
rescatters off the quark through the direct and crossed 
channels (see Fig. This interpretation is possible be- 
cause the outgoing hadron and the hadron appearing in 
the loop are the same. In general, however, it is too 
restrictive to treat final state interactions exclusively as 
a reinteraction of the outgoing hadrons. In fact, in an 
interference-type cut-diagram, any kind of final state in- 
teraction leading to an imaginary part in the amplitude 
of the process q* hX generates a contribution to the 
Collins function. In our case, for instance, it is sufficient 
to employ a different particle in the loop to generate a 
final state interaction which is not a rescattering of the 
outgoing hadron. 

As a final remark, we would like to mention that the 
model we discussed could be applied to describe the 
quark distribution inside a pion in an attempt to gen- 
erate T-odd distribution functions [^0. However, one 
can readily see that, because of the different kinematical 
conditions, self-energy and vertex diagrams do not ac- 
quire an imaginary part, which is essential for producing 
a non-zero T-odd function. 

In conclusion, we have shown that the Collins function 
can be generated using a simple pseudoscalar coupling 
between quarks and pions to model the fragmentation 
process. The Collins function turns out to be zero at the 
tree level due to the absence of any final state interaction 
which is at the origin of its existence, but renders itself at 
the one-loop level. The calculation performed here sug- 
gests that the Collins function, being non-zero already in 
a simple model, is very unlikely to vanish in reality. It is 
therefore a worthwhile task to pursue the measurement 
of in DIS and e+e~ annihilation. Moreover, our re- 
sult supports the idea of using one-particle inclusive DIS 
as a promising process to investigate the transversity dis- 



tribution of the nucleon. 
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